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Within the t-J model, we study the electronic properties of the doped antifer- 
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romagnet on the kagome lattice based on the framework of the self-consistent 
£^ . mean-field theory. At the half-filling, the spin-liquid ground-state energy per 

o, 

site of the kagome antiferromagnet is E g /N s J = —0.859, which is in very 
good agreement with the numerical estimates. Away from the half-filling, the 

electron photoemission spectroscopy and density of states are discussed, and 

O . 

\Q . the results indicate that there is a gap in the normal-state of the system. 
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It has been clear in the last decades that detailed understanding of the motion of the 
holes in a two-dimensional (2D) antiferromagnet is one of the central issues in the theory 
of strongly correlated electron systems, especially as it relates to the copper oxide super- 
conductors 0. This followed from the argument, made by many researchers ||, that the 
essential physics of the copper oxide materials is contained in 2D doped Mott insulators, ob- 
tained by chemically adding charge carriers to a strongly correlated antiferromagnetic (AF) 
insulating state, therefore the physical properties of these systems mainly depend on the ex- 
tent of doping, and the regimes have been classified into the underdoped, optimally doped, 
and overdoped, respectively. The normal-state properties of the copper oxide materials in 
the underdoped and optimally doped regimes exhibit a number of anomalous properties in 
the sense that they do not fit in the conventional Fermi-liquid theory ||. In particular, 
the unusual normal-state gap behavior ||, and the stripe order of holes and spins || are 
found in the underdoped regime. It is also shown from the experiments || that the high- 
est superconducting transition temperature T c for the square lattice Cu0 2 planes in the 
tetragonal structure in the optimally doped regime is reduced by the structure distortions 
of the Cu02 planes in the orthorhombic structure in the underdoped regime, and the stripe 
order of holes and spins is closely related with this distortion of the lattice ||. The spinon 
frustration is also induced due to the distortion of the lattice of the system. Motivated 
by the above experimental facts, we believe that the doped kagome lattice antiferromagnet 
(KLA) is an attractive candidate for describing the distorted square C11O2 planes since it 
has the same coordinate number with the square lattice. On the other hand, the doped 
KLA, which is the system with the geometric frustration as in the triangular lattice, is also 
of the theoretical interest in their own right, with many fascinating theoretical questions 
remaining unanswered. The undoped KLA has been extensively studied f?|-|T0|j, and the 
results show that its ground-state is the novel quantum disordered spin-liquid state without 
the AF long-range-order (AFLRO). It is also believed that the undoped KLA involving only 
nearest-neighbor coupling is the unique model to display the quantum spin liquid ground- 



state pl| . Recognizing that these novel physics showed by the undoped KLA, a natural 



question is what is the electronic properties with doping on this system. 

The situation for the doped square lattice antiferromagnet (SLA) is a bit more ad- 
vanced. Some physical properties of the doped SLA have been studied intensively by many 



researchers within the t-J model [12], where the strong electron correlation manifests itself 



by the electron single occupancy local constraint, and therefore the crucial requirement is 
to impose this local constraint in analytical calculations [[H|. Based on the charge-spin 
separation, it has been shown that a useful method to treat this local constraint is the 
fermion-spin theory Jill. The fermion-spin theory has been used to discuss the physical 



properties of the doped SLA |T5[ , and the results are consistent with the experiments and 



numerical simulations. In this paper, we apply the fermion-spin approach [14[] to study the 
physical properties of the doped KLA. Our results show that at the half-filling the ground- 
state of the kagome Heisenberg antiferromagnet is the quantum spin liquid state with energy 
per site E g /N s J = —0.859, which is in very good agreement with the numerical simulations 
Jr] [Kj . Away from half-filling, the electron photoemission spectroscopy and density of states 
are discussed, and the results indicate that there is a gap in the normal-state of the system. 
The simple way to visualize the kagome lattice is to regard the triangular lattice as 
consisting of four sublattices and remove the spins on one of the sublattices, therefore there 
are three inequivalent sublattices A, B, and C even without AFLRO, where the spins in A, 
B, and C sublattices pointing to the vertices of an equilateral triangle are placed with no 
two nearest-neighbors pointing in the same direction, then the kagome lattice structure is 
much more complicated than the square lattice. For convenience, the equilateral triangle 
ABC is chosen as the cell as shown in Fig. 1, in this case the kagome lattice is reduced as 
the triangular structure, then the position of the spin s in the cell i is specified by the vector 
X is = Ri + d s , with s — A, B, C, i.e., each cell contains three spins, while the lattice vector 
Ri and the reciprocal lattice Kj satisfy the relationship as Ri ■ Kj = 2ir5ij. With the above 
definition, the 2D t-J model on the kagome lattice is described by the Hamiltonian as, 

H = —t 2^ C isCr C is+VsjCr + h.c. — n2_^C isCr Ci 3i(T + J 2_^ Si s ■ S ie+Vs , (1) 

is,Vs,cr i s ,cr i s ,Vs 



where Cl sff (Ci stff ) are the electron creation (annihilation) operators, Si a = C ia aC is /2 are 
spin operators with a = (a x , cr y , a z ) as pauli operators, and \i is the chemical potential. The 
sum is over all sites X is , and for each X is , over the nearest-neighbor r] s . The Hamiltonian 
(H) is restricted to the sub-space where a given site cannot be occupied by more than 
one electron, i.e., X^ Cj s ff Cj aj(T < 1. In the fermion-spin representation [TJJ], C^ = h\S^, 
Cii = hlS^, where the spinless fermion operator hi keeps track of the charge (ho Ion), while 
the pseudospin operator S^ keeps track of the spin (spinon), the t-J model ([[]) is decoupled 
within the mean-field approximation (MFA) as [[HJ H = H t + Hj — 24Ntx4> with 

H t = 2-xt £ hl s+v h is -^hlh is , (2a) 

H J = 2 Jeffe 5Z ( S t S *s+ris + S i. S t+ri.) + J eff J2 S t S L+Vs' ( 2b ) 

where N is the number of cells, J e // = J[(l — 5) 2 — 2 ], and e = 1 + 2t(j)/ J e ff. The nearest- 
neighbor spin bond-order amplitude x an d holon particle-hole parameter <p are defined as 
X = (Si^S~ +rjg ), (j) = (h ig h is+ris ) , respectively. In the present mean-field case, the holon 
moves in the background of spinons, and the spinon's behavior is described by the anisotropic 
Heisenberg model. Since there are three inequivalent sublattices A, B, and C in the kagome 
lattice system, then the one-particle spinon and holon two-time Green's functions are defined 

as matrices, 

fD AA (k,t-t') D AB (k } t-t') D AC (k,t-t')\ 

D x (k,t-t')= D BA (k,t-t') D BB (k,t-f) D BC (k,t-t') , (3) 

\D^(k,t-t') D c x B (k,t-t>) D c x c (k,t-t')) 

where A = s, z, h, L>f (V - j v ,t - t>) = ((S+(t);S-(t'))) and Df(v ~ >»* " *') = 

((Sf (i); Sj u (t'))) with /i, v = A, B, C, are the spinon's propagators, while D^ v {i^ — j u ,t — 

t') = ((hi (t); hj v {t'))) is the holon's propagator. At the half filling, the t-J model is reduced 

as Heisenberg model. Since the absence of the simple three sublattice magnetic order on 

KLA has been convincingly demonstrated by many researchers |TT| , |T6|| , then in the following 

discussions, we will study the holon moves in the disordered spin liquid state, where there is 

no AFLRO, i.e., (Sf ) = 0. In this case, the basic equation for the Green's functions within 



t-J model on the square lattice has been discussed in detail in Ref. |fL5j. Following their 



discussions [jT5|] , the Green's functions (3) is obtained as, 

fT^{k,u)a^(k) T^(k,Lo)aj 2 (k) T^(k,u)a} 3 {k)\ 
T^{k,u)a^{k) T%{k,u)a%{k) T%{k,u)a%{k) 



D x (k,uj) = Y, 

3=1 



X 



Vr&fo «)<£(*) r^(k,u)a^(k) T^(k,uj)aj 3 (k)J 

1 



[^ - ^(k)]^ - u 2 2 (k)][^ - u 2 3 (k)l 
where T^(k,u) is the adjugate matrix of 



( u 2 - (%l + %;7fci) ^37fc4 - Z&Tikh 



A^(k,u) 



U! 



z&iU - z viU u 2 - {z ix + z i2lk2 ) 

V z &l*m - z &lii z ^lts - z £X1m 

with £ = s, z, and T h (k,uj) is the adjugate matrix of 

/ uj + fi -4%t7 fc4 -4xt7 fc6 \ 



^37fc6 - Z&Jk7 \ 
Z^3lk8 — Z^lk<A 



Z& + Zj^rik?) J 



A h (k,u>) 



-^xti/U u + 1* ~^Xtlk8 

V-4xt7fc 6 -4%t7fcs w + n J 



while, 



(4) 



(5) 



(6) 



a?(k)= % 7fc * 4 -2% Z i5lk8 , (7) 

V Z^lla Z^Yks -2% / 
and a h = 1, where Z sl = 8J| // [e 2 (aC + 1/2) + (aC z + 1/2)], Z s2 = 4J 2 // (e 2 ax 2 + eax), 
Z s3 = AJj ff [e(aC + aC z + 1) + 3(e 2 a% + eaXz) - (e 2 «x 2 + ea%)], ^4 = 4J e// (ex + X*), 
Z s5 = 4J e// (e X2 + x), Z zl = 16J 2 //e 2 (aC + l/2), Z z2 = 8J 2 ff ea X , Z z3 = Z 21 /2 + 2Z 22 , Z z4 = 
4:J eff x, Z z5 = Z zi , 7 fcl = cosky+cos(k y -k x ), 7 fc2 = cos/c^+cos/cj,, 7 fe3 = cos^+cos^-A^), 
7fc4 = (l+e ifc »)/2, 7fc5 = (e^+e J fe-^))/2, 7fc6 = (l+e^- fe *))/2, 7fc7 = (e"^+e^)/2, 7fc8 = 
(1 + e~ lkx )/2, 7 fc9 = (e -4 ^ + e i ( kv ~ kx ^)/2, and the order parameters C — J2 (^t+v^i's+v /)> 
Xz = 2(Sf s Sl +rjs ),C Z = 2 J2 (Si s+T]a S i s + Vs , ) ' while U V ( k ) is the solution of the determinant 
\A\(k, u)j)\ = 0. In order not to violate the sum rule (S^S~) = 1/2 in the case (£?) = 0, the 
important decoupling parameter a has been introduced in the above calculation, which can 



be regarded as the vertex corrections [|15| , |l7P . With the help of the spectral representation 
of the correlation functions, the decoupling parameter a and order parameters \i Xzi C, C z , 
<f), and chemical potential ft can be determined by the seven self-consistent equations ||15| . 

At the half-filling, the t-J model is reduced as the Heisenberg model, where e = 1, 
Xz — X-, Cz — C in the rotational symmetrical case, and therefore D z (k,u) = D s (k,uj)/2. In 
this case, the disordered spin liquid ground-state energy per site is E g /N s J = —0.859. For 
comparison, some results of the ground-state energy of KLA obtained from the numerical 
simulations and the present theoretical result are listed in Table I, where we find that the 
present spin liquid energy seems to be closer to the extrapolated finite lattice results of 
Leung and Elser ||, only is more than 3% higher than the best numerical estimates by Zeng 
and Elser , and is almost identical with that of the variational calculation by Sindzinger, 
Lecheminant and Lhuillier |IJ . Therefore the accuracy of the present mean-field results at 



half-filling is confirmed. As a by-product, the Heisenberg antiferromagnet on the triangular 
lattice has been discussed, and the spin liquid energy per site is E g /NJ = —0.966, which 
is essentially identical to results obtained by Kalmeyer and Laughlin |1| and Lee and Feng 
|T9T] based on the resonating-valence-bond state. 

Away from half-filling, we find that the phase separation, present in the doped SLA |20 ], is 
absent here. For discussing the physical properties of the electronic state of the doped KLA, 
we need to calculate the electron Green's function G^(i^ — j u , t — t') = ({C\ CT (t); Cj va (t'))), 
which is a convolution of the spinon Green's function D s (k,t — t') and holon Green's function 
Dh(k, t — t') in the framework of the fermion-spin theory, and can be obtained at the mean- 
field level as, 



1 roo dtu' f 00 du" n F (uj') +n B (u") 



( ' A£»(p-k,u>)A: A (p,J>)^ 



g -(^)^e/:^/_ 



. . 27T UJ + LO 1 - UJ" 
P 



A^(p-k,cu')A» A (p,u;") 
A A »(p-k,uj')A» A (p,u»)) 



where A^ u (k, u) = — 2ImD^ 1/ (/c, uj) and A^ u (k, to) = —2ImD^ u (k, uj) are the holon and spinon 
spectral functions, respectively, n B {uj) is the boson distribution function, while n F {uj) is the 



6 



fermion distribution function. From the Green's function (H), we obtain the electron spectral 
function as A{k,w) = -2ImG AA (k,u) - 2ImG BA (k,iu) - lmG CA (k,u). 

Although the structural distortions of the Cu02 planes in the copper oxide materials 
in the underdoped regime, measured most conveniently by the copper-oxygen bond angles 
H, are very complicated and materials dependent, the common feature is that the spin 
frustration has been introduced. This spin frustration is related to the variation of the 
electronic density of states at the Fermi energy || . Among the doped antiferromagnets the 
most helpful for discussing the effects of the spin frustration on the electronic properties due 
to the distortions of the square C11O2 planes may be the doped KLA since, as mentioned 
above, it has the same coordinate number with the square lattice. We have performed the 
numerical calculation for the electron spectral function A(k, u>) of the doped KLA, and the 
results at the doping 5 = 0.12 for the parameter t/J = 2.5 at the zero temperature is plotted 
in Fig. 2. In comparison with the results of the spectral function on the square lattice ||15|| , 



it is shown that the structure of the spectroscopy on the kagome lattice is much complicated 
than the square lattice. Moving in the momentum space in the direction of O point, the 
photoemission spectroscopy weight increases, while the inverse photoemission spectroscopy 
weight decreases. On the contrary, increasing momentum towards the M point, the situation 
is inversed. For the further understanding of the electronic state properties, we investigate 
the electron density of states p(w), which is closely related with the spectral function, and 

is given as p(uj) = jjJ2A(k,uj). The results of p(uj) with (a) the doping 5 = 0.12 and (b) 

k 

the doping 5 = 0.06 for the parameter t/J = 2.5 at the zero temperature are plotted in 
Fig. 3, where the existence of the gap in the electron density of states is an important 
feature. It is also shown that the density of states is shifted towards smaller energies with 
increasing doped holes, and the total weight is reduced since the integral of the density of 
states up to the Fermi energy has to be equal to the number of electrons, therefore there 
is a tendency that the gap narrows with increasing dopings since some states appear in the 
gap upon dopings. These results are consistent with the experiments of the copper oxide 
materials in the underdoped regime 0]. In comparison with the results of the density of 



states on the square lattice f!5 |, we find that the electron density of states near the Fermi 
energy on the kagome lattice is lower than these on the square lattice. Irrespective of the 
coupling mechanism responsible for the superconductivity in the copper oxide materials, the 
transition temperature T c is usually related to the electronic density of states at the Fermi 
energy. Our results may interpret that the highest superconducting transition temperature 
T c for the square lattice CuO<i planes in the tetragonal structure in the optimal doping is 
reduced by the structure distortions of the Cu02 planes in the orthorhombic structure in the 
underdoped regime ||, i.e., this observed decrease of T c in the underdoped regime || is due 
to the split singularity of the electronic density of states with the gap and lower electronic 
density of states near the Fermi energy. Moreover, the present results also indicate that the 
existence of the normal state gap and decrease of the electronic density of state near Fermi 
energy is possibly induced by the distortions of the lattice (then the spinon frustration). 

In summary, we have discussed the electronic properties of the doped KLA within the t-J 
model. In the self-consistent mean-field level, we show that the ground-state of the undoped 
KLA is the spin liquid state with the energy per site E g /N s J = —0.859, which is consistent 
with the numerical simulations |7|-|iT)[ . It is also indicated that there is the normal-state gap 
in the electronic density of states for the doped KLA. 
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FIGURES 

FIG. 1. The spin configuration of the sublattices A, B, and C on the kagome lattice. 

FIG. 2. Electron spectral function A(k,uf) of the t-J model on the kagome lattice within the 
self-consistent mean-field theory in the doping 5 = 0.12 for the parameter t/J = 2.5 at the tem- 
perature T = 0. Note that there is a gap in the electron spectroscopy. 

FIG. 3. Electron spectral density of the t-J model on the kagome lattice in the doping (a) 
5 = 0.12 and (b) 5 = 0.06 for the parameter t/J = 2.5. 
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TABLES 
TABLE I. A comparison of the ground-state energy per site for the antiferromagnetic Heisen- 

berg model on the two-dimensional kagome lattice. 



Authors 


E g /N s J 


Method 


Zeng and Elser @ 


-0.882 


Finite lattice 


Yang, Warman and Girvin || 


-0.788 


Variational Monte Carlo 


Leung and Elser || 


-0.877 


Finite lattice 


Sindzingre, Lecheminant and Lhuillier [[Hj 


-0.84 


Variational Monte Carlo 


The present work 


-0.859 


Green's function method 
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